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Abstract 

The order a 4 Roo corrections to positronium P levels are found. The calculation 
is reduced to the ordinary perturbation theory for the nonrelativistic Schrodinger 
equation. The perturbation operators have the Breit-type structure and are ob- 
tained by calculating on-mass-shell diagrams. Found energy corrections constitute 
numerically 5E{2 1 P l ) = 0.06 MHz, 6E(2 3 P 2 ) = 0.08 MHz, 5E(2 3 P 1 ) = 0.025 MHz, 
8E(2 3 P 2 ) = -0.58 MHz. 



1 Introduction 



Precision measurements of positronium structure provide a unique test of quantum elec- 
trodynamics. The typical accuracy reached in the measurements of the positronium 2 3 S\ — 
2 3 Pj(J = 0,1,2), l 3 Si - 2 3 Si and 2 3 S 1 - 2 1 P 1 intervals constitutes few MHz Q, @, [| 
|, | |, 0. The order a 3 log(l/a)-Roo and a 3 ^ corrections (i^ = 109 737.315 682 7(48) 
cm -1 is the Rydberg constant) || BL |10| are insufficient now for the comparison of quantum 
electrodynamics with those experimental data. 

The two-body bound state QED problem is by itself of independent theoretical interest. 
The generally accepted theoretical approach to it goes back to Refs Jll|, |l|, This 



approach starts from the introduction of a relativistic two-body wave equation, which can 
be solved exactly and in the nonrelativistic limit reduces to the Schrodinger equation. Then 
a perturbation series is developed about the exact solution. 

Our approach is different and is described, as applied to the corrections logarithmic 
in a, in Refs |TJ], |1|, jlB]] . The corrections discussed in Refs |TJ], [15], [[(J and here, are of 



relativistic origin and can be found as follows. We construct effective perturbation oper- 
ators by calculating on-mass-shell diagrams and then use those operators in the standard 
perturbation theory for the nonrelativistic Schrodinger equation. 

The order a 4 log(l/a)i?oo corrections to the positronium levels were calculated recently 



by Fell [17j and then by us |I| (we have corrected in Ref |L6| a numerical error made in Ref 
fT5|j). This shift exists in the S states only and scales with the principal quantum number 
n as n~ 3 . The logarithmic structure of this correction allows one to treat the relativistic 
effects as a perturbation when deriving the result. However, if one tries to go beyond the 
logarithmic approximation, the logarithmic integrals which are cut off at the electron mass 
m should be treated exactly, and the problem becomes extremely tedious. 

Fortunately, for states of higher angular momenta, L > 0, the situation is better since 
their nonrelativistic wave functions fall off at small distances. Therefore, the integrals 
arising in the perturbation theory, converge in the nonrelativistic region which makes the 
problem quite tractable. The main complication (underestimated by us at the beginning 
of the work) is of the "book-keeping" nature. 

Here we present the results of our analytical calculations for the order a 4 i?oo corrections 
to positronium P levels. In the case n = 2 the results for the fine-splitting of P levels can be 
directly compared with the data extracted from the experimental results of Refs 0, 0, EL [/J. 

Analogous corrections for the electron-electron interaction in helium were obtained nu- 



merically in Ref [IT8I . 



2 Contributions of irreducible operators 



l 



2.1 Relativistic correction to the dispersion law 



Let us start with the kinematic correction generated by the t> 4 /c 4 term in the dispersion 
law for electron and positron, 

2 4 fi 
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By means of motion equations we transform the expectation value to 
(i) ma G I x p 2 1 \ 

H^(s n 3 + ^ + ^ + - + —); S n - 2K - 1 (3) 
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After extracting the overall factor ma 6 everything else here and below is written in the 
usual atomic units. The substitution of the nonrelativistic Coulomb expectation values for 
l/r fc , 

= v ^ " 5?- (4) 
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reduces this energy correction to 
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where e n = ma 6 /n 3 . 



^ = ^-T(8-dr + ^k). ( 6 ) 



2.2 Relativistic corrections to the Coulomb interaction 

This perturbation operator will be extracted from the scattering amplitude for free particles. 
It is convenient to consider the positron as electron of opposite charge. Then the scattering 
amplitude due to the single Coulomb exchange is 

Aixoi 

a c = — jrp(p'> p)p{-p', -p)> ( 7 ) 

where 

p(p',p) =u + (p')u(p), q = p'-p. (8) 

We define scattering amplitudes with the sign opposite to the standard one, so that it 
can be immediately identified with the matrix element of an interaction operator in the 
momentum representation. Let us substitute into (|7j) the solution of the free Dirac equation 



u{p) = J^-\ + (£)w. (9) 
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where w is a bispinor describing particle at rest, 



Kip) 



The corrections of fourth order in v/c are 

Vc = -y^4~5" {(P 2 -P 2 ) 2 ~ ¥,P' x p){q 2 + 3p' 2 + 3p 2 ) - 2(<S,p' x p) 2 } . (10) 

We have neglected here the operator proportional to p 2 + p' 2 since its expectation value 
in the coordinate representation (A<5(r) + 5(f) A) vanishes for P states. The expectation 
value of operator (|H]) is again conveniently calculated in the coordinate representation. Its 
spin-independent part is 

I An 



t(p -v 



„2\2 



4(r- 4 ). 



(11) 
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The treatment of the spin-dependent operators is somewhat more complicated. Here we 
pass over to the Fourier transforms of the operators: 
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then use motion equations and expectation values 
wave function derivative at the origin 
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In this way we get 
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as well as the value of the radial 
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Now, due to the Coulomb interaction electron (positron) can go over into a negative- 
energy intermediate state. The corresponding contributions are described by Z-diagrams 
of the kind presented in Fig.l. Here the particle staying in a positive-energy state can be 
treated in the nonrelativistic approximation. The point is that the large energy denomi- 
nator, equal approximately to 2m ^> E n , due to "heavy" intermediate states, and small 
matrix element of the Z-line provide the sufficient power of a in the perturbation. In this 
way we get for the perturbation operator 

(47ra) 2 f d 3 k 1 



v ( 



c- 



2m J (2tt) 3 k 2 (k - q) 2 

op' ap Bp' a(p + k) a{p + k) op 
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(17) 
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Going over to the coordinate representation, we note that the last integral is in fact the 
convolution of the Fourier-transform of the operator in/r 2 with itself. So, 
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(!) = ma = _ 

c ~ 2 6 \ / 2 5 ■ 5 V 3n 2 1 



Let us note that actually the integral in (|1T|) diverges at large k. It can be easily seen 
however that the divergent part is independent of the momentum transfer q. Therefore 
in the coordinate representation the corresponding operator is just 8{r). Its expectation 
value does not vanish in S states only and leads there (via the accurate cut-off of the linear 
divergence at k ~ m) to a correction of the order ma 5 . 



2.3 Single magnetic exchange 

In the noncovariant perturbation theory the electron-positron scattering amplitude due to 
the exchange by one magnetic quantum is 

4xa j i (p , ,p)j j (-p',-p) ( qiqA . . 

= ~ e, V li - -¥) ■ (19) 

Here 

j(p',P) = u + (p')au(p) (20) 

is the matrix element of the current taken over the solutions @ of the free Dirac equa- 
tion. In the dispersion law for electron and positron it is sufficient here to confine to the 
nonrelativistic approximation. 

Let us start from the contribution to the perturbation operator produced by the v 2 /c 2 
corrections to the currents: 

V ~ a / V 2 (v 2 - + nM^^l + 1 (ffS) 2 

+ ^ ((,' 2 - P 2 ) (l + - *j + (2p + gl S (qS) ] } . (21) 

Here we have neglected again the terms with vanishing P state expectation values. Going 
over to the coordinate representation by means of (|13|), ( fUD and 

d 3 q k-KqiUjjqr = S Jj ~ n i n j ^ 



(2tt) 3 g 4 2r 



,iqr 
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(23) 



we obtain 



(24) 
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Let us consider now the retardation effect. To this end the currents can be taken in the 
leading approximation: 

j(p',p) -> — {$' + P + iq* a). (25) 

At the atomic momentum transfer, q ~ ma, the perturbation of interest originates from 
the second-order term of the expansion of the factor [E n — (p 2 + p' 2 )/2m — q]^ 1 in ( |19D in 
powers of (E n — (p 2 + p' 2 )/2m)/q: 

v -> - -^ (E "V^ )2 {«" + 2 ^' - 2 ™ - 2i ^ x ^ - ^ + <^ 2 } • 

(26) 

At first sight, the expectation value of this operator diverges linearly at small q. This 
divergence can be demonstrated however to be unrelated to the order ma 6 correction we 
are interested in. Indeed, let us split the region of integration over q into two, those from 
to A and from A to oo where ma 2 C A < ma. In the second region our expansion is 
applicable and the result of integration contains a term proportional to 1/A. Since the initial 
integral is independent of A, this term cancels in the sum with the integral over the first 
region calculated without the expansion. Meanwhile, this last integral has no contribution 
of the a A Roo order independent of A. Therefore, taking as the Fourier-transform of 

£ (ar _ Ml) , (27) 

the operator 



-^(Sp 2 -(np) 2 ), (2* 



we get finally: 



-,(1) « i .. , 15 13 



E '« = W^{ U -n + M)- (29) 



Magnetic quantum propagates for a finite time and can cross arbitrary number of the 
Coulomb ones. Simple counting of the momenta powers demonstrates that it is sufficient 
to include the diagrams with one and two Coulomb quanta (dashed lines) crossed by the 
magnetic photon (wavy line). In the first case, Fig.2, the perturbation operator arises as a 
product of the Pauli currents ( |25| ) and the first-order term in the expansion in (E n — (p 2 + 
p' 2 )/2m)/q: 

a 2 f d 3 k (4tt) 2 ( pZ+p'2-gk + k 2 
Vm ° = m 2 lj2^^W\ ^ K]X m 

+ 2pp' - 2 M1M _ 2i(p xk,S)- k 2 f + (kS) 2 
k 2 

In the second case all the elements of diagram 3 should be taken to leading nonrelativistic 
approximation: 

a 3 r d 3 k r d 3 k' (47r) 3 
MCC = ~2m~ 2 l (2^7 (2ny k *(q-k>) 2 (k'-k) 2 X ■ 

k 2 + 2 pp' - 2 MKM _ 2%{p x k, S) - k 2 S 2 + {kSf 
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All the integrals in (|30D and ( j3lD constitute convolutions of the already known Fourier- 
transforms of powers of r. In this way we get for the first operator: 



E^) n = — — - — - \ — 13 H 



30 13 



JMC 2 5 • 3 • 5 I n n 2 

3n 2 J 5 V n 



3SL ( 2 \ (SL) 2 -AS 2 f 2 . . 

' 29 + — ) + 2^^ ( 13 - - ) } ; (32) 



and for the second one: 

E { Scc = i^rtl { 15 ( l - I) + ^ L - 2 (^) 2 + 8 ^} • ( 33 ) 

One more energy correction of the a 4 Poo order at the single magnetic exchange is due 
to Z-type diagrams (see Fig.4). To leading approximation one gets easily 

v »- = ^lwf¥^W { kl ~ ' m ~ f) + ~ kS ' {? ~ x ' (34) 



Standard calculations give now 

_Jn 

2 2 • 5 2 1 2 ' 4~~ v ~~ y ~ ) V 3n 



E { mI = -^\" + iSL-(SL) 2 'S 2 }(i-^- 2 ). (30) 



2.4 True radiative corrections 

Curiously enough, even the true radiative corrections of the a 4 Poo order to the P-levels 
energy can be presented in a simple form practically without special calculation. They can 
be easily demonstrated to be confined here to anomalous magnetic moment contribution to 
the single magnetic exchange, i.e. to a trivial modification of the usual Breit Hamiltonian. 
With this contribution the Pauli current becomes 

J(p',p) -»• —ip' +p + igqx a), (36) 
2m 

where 

a /a x 2 



= 1 + — -c(-j +..., 0.328 (37) 

Let us note here that the anomalous magnetic moment contributions to the first-order 
retardation effect and to diagram 2 cancel. It corresponds to the absence of the order v/c 
corrections to the Breit magnetic exchange. In this way we get the following correction to 
the P-level energy: 



E ra d = -^f- 5 -2{ U + 4C J SL+(l-8c)^-(SLY -S j). (38) 

Two its features are noteworthy. First, it vanishes in singlet states which looks rather 
natural. Second, it has an extra factor l/V 2 which reflects its radiative origin as distinct 
from the relativistic origin of other order a^R^ corrections to P-levels. 
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2.5 Double magnetic exchange 



Let us consider now irreducible diagrams with two magnetic quanta. To our approximation 
they are confined to the type presented in Fig. 5. Their sum reduces to 



a 



d 3 k (4vr) 2 J (kp)(kp') (kp)(kk')(kp') kk' 



Vmm ~ j^M^Y f '~ 2 p + 

+ iS I qX V- k x p- | }> > (39) 

where fc' = q — k. Now we substitute the products of operators in the coordinate rep- 
resentation for the convolutions of their Fourier-transforms and take the corresponding 
expectation values. In this way we get 

<>„ = ^{l0(l-i)-35L(l-A)}. (40 ) 



3 Corrections of second order in the Breit 
Hamiltonian 

Next class of the order a 4 i?oo corrections originates from the iteration of the usual Breit 
Hamiltonian V of second order in v/c. 

Before proceeding to the calculation, let us note that the Breit potential for the positro- 
nium (see, e.g., [HJ, §84), 

V = -^ + ^S(f)-—^(p 2 + (np) 2 ) + —-^LS 
Am 6 m z 2m z r v J Zm^r 6 

(«5n) 2 --6 ^ +-? -2 U(r), (41) 



2m 2 r 3 I 3 J m 2 

has nonvanishing matrix elements for |AL| = 0, 2 only. A contribution of virtual F-levels 



will be found later. And now we average the angular dependence of (41) over a P-state 
wave function: 

,2 -a 3 / \2 4 ->2 N 



3(<Sn) -5 ->-~(2(£l) +5L--5 . (42) 



After this procedure the perturbation ( ^1 ) can be presented as follows: 



ma 4 



V - (43) 

v = (a{h,^} +b[h,ip r ] +c^) , (44) 

where 

a = -3, 6 = 1 + k, c = -4 + k; (45) 

65L - 3(SL) 2 + 2S 2 

K = ; (46) 

5 

7 



p r = —i(d r + 1/r) is the radial momentum, while 

k = f + l-i 

2 r r 

is the unperturbed Hamiltonian for the radial motion with L = 1 in the Coulomb field. 
According to the standard rules, 

£<?> = (VGV) = ^(vgv), (47) 

G = E .TO = ^_,, (48 ) 
^ E n - E k ma 1 

9 = ^mm, (49) 

Since relativistic effects are included here into the perturbation, the intermediate states 
\kP) are merely eigenfunctions of the nonrelativistic Schrodinger equation in the Coulomb 
field. 

The representation ([H]) enables us to find the energy correction fl4"T|) without recourse 
to the exact form of the Coulomb Green's function g. Indeed, 

(vgv) = - ((2a£ n i - a \h, ;] + b[h, ip r ] + c^) gv 

+vg {2aE n \ + a [h, J] + 6[/i, zp r ] + c4r) ) (50) 
= -afi„ (2 («>+a /J <«>-<S /J «»- |({t;,J}) 

+^<[ip rj t;]) + c(a r <v)-(a r «)). (51) 

When passing from (|50"D to ([51]) we used the equation of motion and presented the pertur- 
bations 1/r and 1/r 2 as the result of (3 and 7 variations in the modified Hamiltonian, 

I r A r 

Derivatives in fl51"|) are taken at /3 = 7 = 1. 

Substituting flHj) into (|5l| ) provides us with the expression containing just the mean 
values of l/r fc (see (§)). 

We collect now coefficients at the powers of 1/n to obtain: 

3a 2 + Uab + 13b 2 2c(2c + 9a + 96) 



(2) = ma" 

p 128n 3 1 15 27 



2c 2 2 /lla 2 + 13a6 + 6o 2 \ 5a 2 , 

-*r + s? 5 + H-2^^- (53) 



Finally, substituting the expressions for a, 6, c (f45|), we get the contribution of intermediate 
P states to the iteration of the Breit Hamiltonian: 

( 2) _ / 1022 - 844k + 227k 2 (4 - k) 2 102 - 29k + 2k 2 45 1 

p ~M 2^ • 3 3 • 5 26-3-n + 26 ■ 5 • n 2 2^3 [ " ( } 



Transitions to intermediate F states are driven by the operator (<Sn) 2 contained in fl4ll) . 
Due to the conservation of the total angular momentum such transitions are possible in 
the states with J = 2 only. The matrix element squared of the operator (Sn) 2 is easily 
calculated by means of closure: 



m\(sft) 2 \F 2 ) 



{P 2 \{Snf\P 2 ) - (P 2 \{Sn) 2 \P 2 f 
(P 2 \(Sn) 2 \P 2 ) - (P 2 \(Sn) 2 \P 2 )' 
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We have used here the identity (Sn) 3 \P) = Sn\P). 

To calculate the radial part of the correction, proportional to 



1 \kF)(kF\ 1 

3 2^ 



8 k r 3 



we present the operator 1/r 3 in the form 



1 

18 
1 

18 



h 

hp 



-iPr 



-ip r 



5 1\ ■ 
r + 5) hF . 

ip r + - - 7 ) — ( -iPr H -) h 



where 



V 
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1 



h F = q- + - - - 



(55) 



(56) 



(57) 
(58) 

2 ' r 2 r (59) 
is the unperturbed Hamiltonian for the radial motion with L = 3 in the Coulomb field. 
The energy correction is again expressed in terms of (r~ k ): 



E 



(2) 



ma 



3/2 7 

2 8 ■ 5 2 XoT 3 _ ^ 
7 

10 --r 



2 6 -5 3 



(60) 



4 Numerical results 

Let us summarize the numerical values of the a 4 Roo corrections to the energies of positro- 
nium 2P levels. They are 

0.06 MHz for 2 1 P l , 

0.08 MHz for 2 3 P 2 , 

0.025 MHz for 2 3 P 1; 

- 0.58 MHz for 2 3 P - 

We wish to emphasize here that the last correction is quite comparable in magnitude 
to the corresponding logarithmic correction (of order a 4 log a R^) to the positronium 2S- 



levels. The latter constitutes for instance 0.96 MHz for the 2 Si state [17, [Lq]. Therefore, 
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there is no special reasons to expect that the nonlogarithmic corrections (of order a 4 i?oo) 
to the positronium 5 levels are small. It makes their calculation quite an actual problem. 

Including a 2 -Roo and a 3 -Roo terms we obtain the fine-structure intervals between 2 2 +1 Pj 
levels. In the Table 1 these theoretical values are compared with the transition frequencies 
extracted from the results of two recent experiments |2|, |], H, [7| (all systematic and statistical 
errors are added quadratically when extracting the experimental numbers). 
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Table 1: Fine-structure intervals between 2 2S+1 Pj levels (in MHz). 





Experiments 


Theory 




Michigan [|. | 


Mainz §, 


Total 


a^-Roo contribution 


E( 3 P 2 )-E(*P ) 
E( a P 1 )-E(^Po) 
EfPj-EfPo) 


9884.5±10.5 
5502.8±10.9 
7323.1±16.5 


9875.27 ±4.44 
5487.23 ±4.50 
7319.65 ±7.64 


9871.54 
5485.84 
7312.88 


0.66 
0.60 
0.64 
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Figure captions 

Fig.l. Z-type double- Coulomb exchange 

Fig. 2. Single-magnetic-single-Coulomb exchange 

Fig. 3. Single-magnetic-double-Coulomb exchange 

Fig. 4. Z-type single-magnetic-single-Coulomb exchange 

Fig. 5. Double- magnetic exchange 
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Fig.l 




